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Abstract 

Prom the current ATLAS and CMS results on Higgs boson mass and decay rates, the NMSSM 
is obviously better than the MSSM. To explain the fine-tuning problems such as gauge hiearchy 
problem and strong CP problem in the SM, we point out that supersymmetry does not need to 
provide a dark matter candidate, i.e., -R-parity can be violated. Thus, we consider three kinds of 
the NMSSM scenarios: in Scenarios I and II i?-parity is conserved and the lightest neutralino relic 
density is respectively around and smaller than the observed value, while in Scenario III .R-parity 
is violated. To fit all the experimental data, we consider the x 2 analyses, and find that the Higgs 
boson mass and decay rates can be explained very well in these Scenarios. Considering the small 
X 2 values and fine-tuning around 2-3.7% (or 1-2%), we obtain the viable parameter space with light 
(or relatively heavy) supersymmetric particle spectra only in Scenario III (or in Scenarios I and 
II). Because the singlino, Higgsinos, and light stop are relatively light in general, we can relax the 
LHC supersymmetry search constraints but the XPNON100 experiment gives a strong constraint in 
Scenarios I and II. In all the viable parameter space, the anomalous magnetic moment of the muon 
(gn — 2)/2 are generically small. With .R-parity violation, we can increase (<? M — 2)/2, and avoid 
the contraints from the LHC supersymmetry searches and XENON100 experiment. Therefore, 
Scenario III with imparity violation is more natural and realistic than Scenarios I and II. 

PACS numbers: ll.10.Kk, 11.25.Mj, ll.25.-w, 12.60.Jv 
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I. INTRODUCTION 



The Higgs boson mass in the Standard Model (SM) is not stable against qunatum correc- 
tions and its square has quadratic divergences. Because the reduced Planck scale is about 
16 order larger than the electroweak (EW) scale, there exists huge fine-tuning around 10 -32 
to obtain the EW-scale Higgs boson mass. Supersymmetry is a symmetry between the 
bosonic and fermionic states, and it naturally solves this problem due to the cancellations 
between the bosonic and fermionic quantum corrections. In Minimal Supersymmetric SM 
(MSSM), the gauge couplings for SU(3)c, SU{2) L and U(l)y gauge symmetries are unified 
at about 2 x 10 16 GeV [T], which strongly suggests Grand Unified Theories (GUTs). Unlike 
the SM, we can have the renormalizable superpotential terms that violate the baryon and 
lepton numbers, and then there may exist proton decay problem. To solve such problem, 
we usually introduce the .R-parity under which the SM particles are even while the extra 
supersymmetric particles (sparticles) are odd. Thus, the lightest supersymmetric particle 
(LSP) like neutralino can be cold dark matter candidate [21 E]- 

However, there are strong constrains on the supersymmetry viable parameter space from 
the recent LHC supersymmetry searches [IHH]. For example, in the Minimal Supergravity 
(mSUGRA) model or Constrained MSSM (CMSSM), gluino mass should be larger than 
about 1.4 TeV and 850 GeV for squark masses around and much larger than gluino mass, 
respectively. Also, squarks (at least the first two generation squarks) must have masses 
larger than about 1.1 TeV from the ATLAS and CMS Collaborations at the LHC @HH] - 

Recently, the ATLAS and CMS Collaborations have announced the discovery of a Higgs- 
like boson with mass around 126.5 GeV and 125.3 ± 0.6 GeV, respectively [9T4TT]. In the 
77 final state, the ATLAS and CMS rates are roughly 1.9 ± 0.5 and 1.56 ± 0.43 times 
the SM prediction. In the ZZ — > 4£ channel, the ATLAS and CMS signals are roughly 
l.lioi ari d 0.7jTal times the SM prediction, respectively. In the bb, t + t~ and WW — > Evii; 
channels, the ATLAS rates are respectively 0.48+2J2, 0.16 jl];;^, and 0.52^^1 times the SM 
prediction, and the CMS rates are respectively 015^, -°- 14 -a73> and q - Q2 -oaI- So these 
rates are somewhat suppressed compare to the SM prediction but error bars are relatively 
large. The Higgs physics implications in the supersymmetric SMs (SSMs) have been studied 
extensively [i~2"HT6] . By the way, the new results from the CDF and DO experiments [17] 
support the ~ 125 GeV Higgs signal and suggest an enhancement relative to the SM of the 



3 



jy+Higgs with Higgs— > bb rate by a factor of 1.97_ ' 68 . But we will consider not it here 
since it is different from the ATLAS and CMS results. 

As we know, there are two Higgs doublets H u and Hd in the MSSM that gives masses to 
the up-type quarks and down-type quarks/charged leptons, respectively. The lightest CP- 
even Higgs boson mass, which is a linear combination of H® and H® and usually SM-like, 
is smaller than Z boson mass Mz at tree level. Thus, to realize the lightest CP-even Higgs 
boson mass around 125.5 GeV radiatively, the squark and/or gluino masses will be about 
a few TeV in general in the mSUGRA/CMSSM. And then there exists at least less than 
one-percent fine-tuning. Moreover, it is difficult to explain the above Higgs decay rates and 
generate the correct Higgs boson mass simultaneously in the MSSM. For example, if the 
SM-like Higgs particle has dominant component from H®, we can suppress the rates in the 
bb and r + r~ final states, and then increase the 77 rate. But the rates for the ZZ — > 4£ and 
WW — > Iviv channels will increase as well. Also, if the stop is light, we can increase the 
Higgs to two photon rate, but it is difficult to generate the 125.5 GeV Higgs boson mass [T2T 
[16]. The possible model might be the light stau scenario [13]. Therefore, we shall consider 
the next to the MSSM (NMSSM) where an SM singlet field S is introduced. The points are 
the following: (1) We can increase the Higgs quartic coupling from the superpotential term 
\SH d H u if the ratio tan/3 = (H®) / (H®) of the vacuum expectation values (VEVs) for H® 
and is not large; (2) We can suppress the couplings between the W/Z gauge bosons and 
the Higgs particle due to the mixings among S, H®, and H®. 

On the other hand, the strong CP problem is another big fine-tuning problem in the SM. 
From the experimental bound on the neutron electric dipole moment (EDM), the strong 
CP phase 9 is required to be smaller than 1CT 10 . An elegant and popular solution to 
the strong CP problem is provided by the Peccei-Quinn (PQ) mechanism [18], in which a 
global axial symmetry U(1)pq is introduced and broken spontaneously at some high energy 
scale. The axion a is a pseudo-Goldstone boson from the spontaneous U(1)pq symmetry 
breaking, with a decay constant f a . The original Weinberg-Wilczek axion [19J is excluded 
by experiment, in particular by the non-observation of the rare decay K — > 71 + a [20J. 
There are two viable "invisible" axion models in which the experimental bounds can be 
evaded: the Kim-Shifman-Vainshtein-Zakharov (KSVZ) axion model [21J and the Dine- 
Fischler-Srednicki-Zhitnitskii (DFSZ) axion model [22]. From laboratory, astrophysics, and 
cosmological constraints, the U(1)pq symmetry breaking scale f a is constrained to the range 
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10 10 GeV < f a < 10 12 GeV [20]. Interestingly, for such f a range, the invisible axion can be 
a good cold dark matter candidate with correct relic density [20J. 

Because axion can be the correct dark matter candidate, supersymmetry may only need 
to solve the gauge hiearchy problem and realize gauge coupling unification. Therefore, we 
consider three kinds of the NMSSM scenarios: in Scenario I, i?-parity is conserved and 
the LSP neutralino relic density is around the observed value; in Scenario II, .R-parity is 
conserved and the LSP neutralino relic density is smaller than the observed value; in Scenario 
III, .R-parity is violated and then the LSP neutralino is not stable. In particular, Scenario III 
is very interesting since it can not only avoid the constraints from the LHC supersymmetry 
searches and XENON100 experiment [23J, but also may relax the other phenomenological 
constraints. Moreover, the proton decay problem can be solved by requiring the baryon or 
lepton number conservation [21], or by requiring the minimal flavour violation [25J. 

In this paper, we shall study the natural and realistic NMSSM. We first briefly review the 
naturalness condition in the SSMs and discuss the NMSSM with and without i?-parity. To 
satisfy the phenomenological constraints and fit the experimental data, we consider the \ 2 
analyses for all three kinds of Scenarios, and find that we can indeed explain the Higgs boson 
mass and decay rates very well. Considering the small \ 2 values and fine-tuning around 2- 
3.7%, we obtain the viable parameter space with light (e.g. less than around 900 GeV) 
supersymmetric particle spectra only in Scenario I. For the small y 2 values and fine-tuning 
around 1-2%, we get the viable parameter space with relatively heavy (e.g. less than about 
1.2 TeV) supersymmetric particle spectra. In particular, the best benchmark point has 
almost minimal x 2 an d 3.7% fine-tuning in Scenario III. The generic features for the viable 
parameter space with smaller x 2 are that the light stop is around 500 GeV or smaller, the 
singlino and Higgsino are light chargino and neutralinos, the Wino-like chargino is heavy, 
and the Bino-like and Wino-like neutralinos are the second heaviest neutralino and heaviest 
neutralinos, respectively. Thus, we find that the LHC supersymmetry search constraints can 
be relaxed due to quite a few jets and/or leptons in the final states in Scenarios I and II, but 
the XENON100 experiment still gives strong constraint on the dark matter direct detection 
cross sections. Moreover, the correct dark matter relic density can be realized in Scenario 
I as well. In particular, tan /3 is not large and the second lightest CP-even Higgs particle 
is SM-like [HJ [TU], which is helpful to increase the SM-like Higgs boson mass. However, 
the additional contributions to the anomalous magnetic moment of the muon (g^ — 2)/2 
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are smaller than three sigma low bound [26] in general due to relatively small tan/3. As we 
know, with i?-parity violation, we can escape the constraints from the LHC supersymmetry 
searches and XENON100 experiment, and the i?-parity violation superpotential term(s) 
may increase the muon (g^ — 2)/2 and explain the neutrino masses and mixings. Therefore, 
Scenario III with i?-parity violation is more natural and realistic than Scenarios I and II. 

This paper is organized as follows. We explain the naturalness criteria in the SSMs in 
Section II. We present the NMSSM with and without i?-parity in Section III, and the exper- 
imental constraints/data and numerical analyses in Section IV. Section V is our conclusion. 



II. NATURALNESS CRITERIA IN THE SSMS 



For the GUTs with gravity mediated supersymmetry breaking, the usual quantitative 
measure Aft for fine-tuning is the maximum of the logarithmic derivative of Mz with 
respect to all the fundamental parameters ctj at the GUT scale 

d\n(M z 



A FT = Max{Af UT } , A GUT 



(1) 



<91n(af UT 

In the following numerical calculations, we will use this definition to calculate the fine-tuning. 

However, the above fine-tuning definition is a little bit abstract. Thus, we shall present 
the concrete bounds on the fi term, third-generation squark masses and gluino mass in the 
following [281 129]. The SM Higgs-like particle h in the MSSM is a linear combintation of H® 
and H®. To simplify the discussion on naturalness, we can reduce the Higgs potential to 

V = ml\h\" + ^\h\\ (2) 

where m\ is negative. Minimizing the Higgs potential, we get the physical SM-like Higgs 
boson mass m^ 

m 2 h = -2m 2 h ■ (3) 
So the fine-tuning measure can also be defined as [28] 

A FT = — ^ • 4 

m 2 h 

For a moderately large tan/3 = (H®)/(H®), for instance, tan/3 > 2, we have 

m 2 h ~ \fi\ 2 + m 2 Hu | tree + m% | rad , (5) 



where /i is the supersymmetric bilinear mass between H u and Hd-, and m 2 H | tr ee and m\ | ra( j 
are the tree-level and radiative contributions to the soft supersymmetry-breaking mass 
squared for H u . Therefore, we obtain the following concrete bounds 



The upper bound on the \i term is 



m h \ ( A FT 



1\ -1/2 



Thus, the /x term should be small than about 400 GeV for 5% fine-tuning. Conse- 
quncely, the charged and neutral Higgsinos will be light. In the NMSSM, we just 
change the \i term to the effective /i term /x e ff = A(S'). 

The one-loop radiative corrections to rn 2 Hu in the leading logarithmic approximation 
from the stop sector are 

3 



5m 



2 I 

ifj stop g7r 2 



!/?K+m| s + IA| 2 )ln(4) , (7) 

where yt is top Yukawa coupling, m~ and m~ are supersymmetry breaking soft 
masses for the third generation quark doublet and right-handed stop, A t is the top 
trilinear soft term, and A is the effective supersymmetry breaking mediation scale. 
Thus, one obtains 

-1/2 , _ , /a -i\ - 1 / 2 



7^ < ! 2 TeV Sin/3 (hi^m ' ( ) 

? 1 + m t 2 ~ u (l+# V 3 J U25.5GevHs%/ 

where x t = At/ \J m \ + m |, > an d £i and are two stop mass eigenstates. Therefore, 



we obtain yrn 2 + < 1.2 TeV. Also, we can require that the lighter sbottom 
mass be smaller than m t ~ 2 , which is automatically satisfied via an simple mathematical 
proof. 

The two-loop radiative corrections to m 2 H in the leading logarithmic approximation 
from gluino are 

*»i.U*- = "J* 2 (v) Win 2 (^) , (9) 
where a s is the strong coupling, and M3 is the gluino mass. Here, the contributions 
from the mixed A t M^ term , which are relevant for large A-term, are neglected. Thus, 
the bound on gluino mass is 



ln(A/TeV)^ ( m h \ ( A FT 



M 3 < 1.8 TeV sin, ( — - ) { ^^) (10) 



1\ -1/2 



So the gluino mass is lighter than about 1.8 TeV. 



Therefore, the natural MSSM and NMSSM should have relatively smaller (effective) \x 
term, stop masses as well as gluino mass. In this paper, we shall not only use Eq. ([I]) to 
calculate the numerical values of the fine-tuning, but also consider the following natural 
supersymmetry conditions: 

• The n term or effective \x term is smaller than 300 GeV. 



• The squar root Mi = \ rri~ + m~ of the sum of the two stop mass squares is smaller 



than 1.2 TeV. Consequencely, we can show that the light sbottom mass is smaller than 



• The gluino mass is lighter than 1.5 TeV. 

However, such kind of the natural MSSM and NMSSM might be excluded by the LHC 
supersymmetry searches and XENON100 dark matter direct detection. Thus, the i?-parity 
violation might be needed for the natural MSSM and NMSSM, and then supersymmetry 
only needs to solve the fine-tuning problem and explain the gauge coupling unification. 

III. THE NMSSM WITH AND WITHOUT /2-PARITY 

Let us explain the convention first. We denote the quark doublets, right-handed up-type 
quarks, right-handed down-type quarks, lepton doublets, and right-handed leptons as Qi, 
D°, Li, and Ef, respectively. We denote the SU(3)c, SU(2)l, and U(l)y gauginos as 
G a , W a , and B, respectively. To solve the \i problem in the MSSM, we introduce a SM 
singlet field S and consider the NMSSM with Z% symmetry which forbids the \x term. The 
superpotential in the NMSSM is 



where yfj, yfj, y\j, A, and k are Yukawa couplings. The effective /i term is obtained after S 
obtains a VEV, i.e., /i e fr = 




Wnmssm 



y%QiU°H u + yfjQiD'H, + y\ 3 UE]E d + \SE d E u + -nS 3 , 



(11) 
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The supersymmetry breaking soft terms are 

3 8 

M 1 BB+M 2 W a W a +M 3 & a G a + H.C. 



a=l 

+m 2 s \S\ 2 + m\\Qtf + m% t |^ c | 2 + m\ c \ D^ 2 + m\\U\ 2 + m\ c \E, 



a=l 

,2 \ttc\2 , 2 



+ m 2 H jH u \ 2 + m 2 Hd \H d \ 2 

2 |7 |2 , _2 I ?c|2 



+ 



Hi' y %l 



y^QiUfHu + yf^Q^H, + yl^L^H, + \A x SH d H, 



i A i 



+-kA k S 3 + H.C. 



(12) 



Similar to the MSSM, the Higgs sector of the NMSSM is described by the following six 
parameters 



A , k , A x , A K , tan/3 , /i efT • 



(13) 



And the supersymmetry breaking soft mass terms for the Higgs bosons rn 2 Hu , rn 2 Hd and m| 
are determined implicitely by Mz, tan/3 and /i e fr via the Higgs potential minimization. 

In addition, from the theoretical point of view, we usually have the family univer- 
sal squark and slepton soft masses in the string model building. Therefore, as in the 
mSUGRA/CMSSM, we consider the following universal supersymmetry breaking soft terms 



.I/, M 2 M 3 - M 



1/2 ? 



(14) 



m 2 ~ 



m| - -2 _ ^2 



(15) 



AU _ Ad 



A% = A . 



(16) 



We consider the NUH-NMSSM in this paper: the Higgs soft mass terms m 2 Hu , m< H d an d m | 
are allowed to be different from Mq (and determined implicitely as mentioned above), and 
the trilinear couplings A\, A K and Aq are not universal. Therefore, the complete parameter 
space is characterized by 



A , k , tan /3 , // eff , A x , A K , A , M_ 



1/2 , M , 



(17) 



where the last five parameters are taken at the GUT scale. 

Next, we consider the i?-parity violation. The most general renormalizable, gauge and 
Z 3 invariant, and .R-parity odd superpotential terms in the NMSSM are 



H^rpv = \SLiH u + -Xij k LiLjEl + \\- k LiQjD c k + -X" ijk U[D c j Dl , (18) 



where Aj, Xijk, A^- fc , and \"ijk are Yukawa couplings. In the above Eq. (18), the first three 
terms conserve the baryon number while violate the lepton number, and the last term 
conserves the lepton number while violates the baryon number. Thus, to forbid the proton 
decay, we require either baryon number conservation or lepton number conservation, i.e., 
we turn on either the first three terms or the last term in the above superpotential [24]. 
The alternative ways are to consider the minimal flavour violation [25J or discrete R- 
symmetry [30]. In particular, the Ay* and A' ijfc terms can contribute to the anomalous 
magnetic moment of the muon (g M — 2)/2 and generate the neutrino masse and mixings, 
and the A^ fe and X'ijk terms can contribute to the b — > 37, etc [21]. We would like to point 
out that the NMSSM with i?-parity violation has been studied before [3T] . and the NMSSM 
with baryon number conservation is similar to the /zi/SSM [32] . 

In this paper, to increase the SM-like Higgs boson mass while keep the sparticle spectrum 
light, we will concentrate on the natural and realistic NMSSM with the following properties: 
(1) tan/3 is not large so that the SM-like Higgs boson mass can be lifted via the tree-level 
\SH d H u term; (2) The second lightest CP-even Higgs boson is the SM-like Higgs particle, 
and then the SM-like Higgs boson mass can be lifted via the mass matrix diagonalization 
from Linear Algebra. However, in such kind of viable parameter space, the muon (g M — 2)/2 
is generically small due to not large tan/3. Thus, to increase muon (g M — 2)/2, we need to 
introduce .R-parity violation \j k and A' ijfc terms in Eq. (18), which will be studied elsewhere. 
Interestingly, we may explain the neutrino masses and mixings simultaneously. 



IV. x 2 ANALYSES FOR THE PHENOMENOLOGICAL CONSTRAINTS AND 
EXPERIMENTAL DATA 

We will consider the x 2 analyses for the phenomenological constraints and experimental 
data in all three scenarios. For our numerical calculations, we use the NMSSMTools version 
3.2.0 [33]. 

In the original package, the points are survived if they satisfy several phenomenological 
and theoretical constraints. Two standard deviation (95% C.L. upper) limits are applied 
for those constraints which have corresponding experimental measurements. In this paper, 
these two standard deviation (95% C.L. upper) limits are replaced by their central values 
and the experimental errors, which are used to construct the global \ 2 ■ There are two 
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advantages for this global fit: (1) The best-fitted benchmark points with minimal x 2 value 
can be found exactly, while the previous method within two standard deviation limits can 
only provide a viable parameter space. (2) The derivation for the central values in the 
two standard deviation limits can accumulate to be a relatively significant drift while the 
global x 2 can have explicit statistical meanings for the 1 or 2 standard deviations from the 
best-fitted points. 

In our analyses, several phenomenological and theoretical constraints are considered. 
These constraints can be divided into the following categories: 

1. The theoretical constraints and phenomenological constraints, which only have 95% 
CL. upper limits, are unchanged in the NMSSMTools 1 . 

2. The following LHC Higgs constraints are added: the second CP-even neutral Higgs 
field i?2 are taken as the SM-like Higgs boson discovered at the LHC and its mass is 
required to be M H . 2 e [124, 127] GeV. All the 5 neutral Higgs fields H u H 2 , H 3 , A 1 and 
A 2 should satisfy the LHC constraints, which are taken as the 95% C.L. on er/ergM 
among the LHC measured Higgs mass regions. Table [I] shows the ATLAS and CMS 
Higgs decay channels that we will consider. 

TABLE I: The LHC collider constraints at 95% C.L. on f/cgM- 7, 8 and 7&8 means LHC 
center mass energy y/S = 7, 8 TeV as well as 7 and 8 TeV combined results. The blank means 
no constraints in this channel. Same conventions are applied in the following tables. "(VH)" in 
the table indicates the experimental results are actually measured in the vector boson associated 
production, which is invariant for the WH production channel in the NMSSM. 

1 For interested readers, the detailed information can be found in NMSSMTOOLS 3.2.0, which corresponds 
to PROB(l)~PROB(29), PROB(35), PROB(41)~ PROB(45), PROB(51) and PROB(52). 
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Channels ATLAS CMS 



H^tt 7, 7&8 

WH -»■ bb 7, 7&8(VH) 

H bb 

H -> ZZ 7, 8, 7&8 7, 7&8 
H -> 7, 7&8 

ff -»• 77 7, 8, 7&8 7, 7&8 
2jfT -> 2j 77 



3. The NMSSMTools two standard deviation constraints are replaced by global \ 2 fits, 
which include: b — > sj, 5m s , 5md, b — > ru T , (g^ — 2)/2 and Br (B — > X s ji + jj,^). We 
update the anomalous magnetic moment of the muon (g^ — 2)/2: Aa M = a° xp — a^ M = 
(28.7 ±8.0) x KT 10 [26]. 

4. LHC Higgs signal strength is constructed in the \ 2 i as shown in Table [h| and Fig. [Tj 
Theoretical predictions in the \ 2 correspond to H2 in the NMSSM. For the 7 and 
8 TeV combined results, the theoretical predicted signal strength for inclusive Higgs 
production channels are combined from the 7 and 8 TeV individual signal strength 
proportional to their accumulated luminosities. 

TABLE II: LHC Higgs signal strength in \ 2 ■ 



Channels 


ATLAS 


CMS 


H —7- TT 


7 


7&8 


WH -»• bb 


7 


7&8(VH) 


H^bb 


7 


7&8 


H -> ZZ 


7&8 


7&8 


H -> W+W- 


7 


7&8 


il — )• 77 


7&8 


7&8 


2jF -> 2.777 


7&8 


7&8 



5. The cold dark matter relic density is 0.112 ± 0.0056 from the seven-year WMAP 
measurements [31]. For the dark matter, we treat it in three different scenarios, as 

12 
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FIG. 1: LHC Higgs signal strength in different production and decay modes. 



shown in the Table III In Scenario I, the lightest stable neutralino is required to have 
the correct dark matter relic density. This is considered in the global x 2 ■ m Scenario II, 
the relic density is required to be smaller than the 95% C.L. experimental upper limit, 
which assumes multi-component dark matter. In Scenario III, the relic density is set to 
be free, which corresponds to the i?-parity violation case. Constraints of effective Higgs 
self-couplings in MicrOMEGAs and spin independent dark matter direct detection of 
XENON100 experiment [23] in the NMSSMTools package are adjusted according to 



different relic density treatment as shown in Table III 



TABLE III: Three different dark matter relic density scenarios. 

scenario I scenario II scenario III 
0,h 2 in x 2 < 95%upper limit free 

Effective Higgs Self-Coupling in MicrOMEGAs > 1 / / x 

XENON100 / / x 
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To be clear, we have considered three scenarios according to the different dark matter 
treatments as explained in category 5. Categories 1 to 4 are common constraints that are 
applied to all three Scenarios. 

The x 2 is constructed as: 

..the _ .exp 

x 2 = B 5i l ) 2 - ( 19 ) 

i 

in which n\ he are theoretical predicted values and fi^ xp are corresponding experimental mea- 
surements. 5i are one standard fluctuations which includes both statistical and systematical 
errors and are taken as the average values for asymmetric errors. 

By adopting the above \ 2 constructions, the number of independent variables in the \ 2 
are: 7 in category 3 (Br (B — > X s {i + fi~) are considered in both low and high dilepton 
energy regions), 14 LHC Higgs decay signal strength (we assume that ATLAS and CMS 
measurements on the same Higgs decay channels are independent), and the different dark 
matter relic density Scenarios in category 5. Besides, there are 9 NMSSM input parameters 



as shown in Eq. (17). So the number of degree of freedom is 22-9=13 for Scenario I and 
12 for Scenarios II and III. The goodness of fit can be shown by comparing the minimum 
X 2 with the n d . 

Note that the current top quark mass m t is 173. 5±1 GeV, we shall choose the central value 
m t = 173.5 GeV in numerical calculations. We emphasize that the SM-like Higgs boson H2 
mass will increase and decrease about 1 GeV if we choose the upper limit m t = 174.5 GeV 
and low limit m t = 172.5 GeV, respectively. Thus, the SM-like Higgs boson H 2 mass range 
from 124 GeV or 127 GeV is fine. Moreover, we define 

R xx = °{PP -> gj) BR(Hj -» XX) 

o{w /ism) BR(/i S m -^XX)' 

where XX can be 77, Z°Z°, W + W~ , bb, and rf. 

We present the B% v versus FCf , RY versus R^ b , Mi versus Mg, A FT versus x 2 , tan (3 versus 
Aa^, versus /x e g, M Hl versus M H2 , Mq versus M g , and Mj versus M fl in Figs. [2j |3j and|4] 
for Scenarios I, II, and III, respectively. The red stars show the best-fitted benchmark points 
with minimal Xmin = 21.16, 19.35, 19.67 for Scenarios I, II, III, respectively. The magenta 



region corresponds to i? 77 > 1.4, Ryv < 1-1, Rbb < 1-0, R TT < 1-0, Mi = . m~ + m~ < 

y ^1 ^2 

1.2 TeV, fi eS < 300 GeV, Mg < 1.2 TeV, \ 2 < Xmin + ^ and A ft < 50. In particular, 
the small \ 2 and Aft < 50 are not compatible with each other in Scenarios I and II due 
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to XENON100 experimental constraint, and then only Scenario III has magenta region. In 
addition to the minimal \ 2 points, we consider three kinds of other benchmark points: 

• Benchmark points IA, IIA, and IIIA have small x 2 , relatively small Aft, the light stop 
mass around 200 GeV, and the first two generation squark masses lighter than about 
1.1 TeV. 

• Benchmark points IB, IIB, and IIIB have small x 2 , relatively small Aft, the light stop 
mass around 200 GeV, and the first two generation squark masses heavier than about 
1.1 TeV. 

• Benchmark points IC, IIC, and IIIC have small x 2 , relatively small Aft, and relatively 
heavier light stop. 

We present the minimal x 2 point, and three other benchmark points in Tables IV, Vj and 



VT| for Scenarios I, II, and III, respectively. Moreover, we study the constraints from the 
LHC supersymmetry searches, and find that only the benchmark points IA and IIA are 
excluded by the current LHC supersymmetry searches [3HHJ . Because the benchmark points 
IXmini IC, and IIC have relatively light spectra (e.g. less thank about 1.2 TeV), the LHC 
supersymmetry search constraints are relaxed in our models, which will be explained briefly 
in the following. The benchmark points in Scenarios I and II have fine-tuning from about 1% 
to 2%, while the benchmark points in Scenario III have fine-tuning from about 2% to 3.7%. In 
particular, in the best benchmark point IIIA we have \ 2 = 21.31, and Aft = 27.0, i.e., 3.7% 
fine-tuning. Also, all the supersymmetric particles are lighter than 830 GeV. By the way, 
all the benchmark points except llXmin satisfy the naturalnes conditions: fi c s < 300 GeV, 
Mi < 1.2 TeV, and Mg < 1.5 TeV. Note that A FT = 130.4 in benchmark point IIxLn; the 
two fine-tuning definitions in Section II are compatible. 

From the viable parameter space in Figs. [2j |3j and [4], we find that tan (3 is generically 
smaller than about 4.5, and then we have the small anomalous magnetic moment of the 
muon (g M — 2)/2, i.e., Aa^ < 4.0 x 10 -10 . Also, we notice the correlation between R^ 1 and 
R\ v which roughly is Rj 1 ~ 1.27 x R\ v . Interestingly, we do have some viable parameter 
space which indeed have R^ 1 > 1.4 and R\ v < 1.1. The generic features for the parameter 
space with small x 2 are that the light stop is around 500 GeV or smaller, the singlino and 
Higgsino are light neutralinos and chargino, the Wino-like chargino is heavy, and the Bino- 
like and Wino-like neutralinos are the second heaviest neutralino and the heaviest neutralino, 
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respectively. Thus, we can understand why the LHC supersymmetry search constraints are 
relaxed: the branch ratios of the first two generation squarks decaying directly to the LSP 
neutralino and quarks are very small around 1%, and the dominant decay channels are 
Wino-like chargino/neutralino and quarks. And then the Wino-like chargino and neutralino 
will decay into quite a few jets or leptons via the light chargino and neutralinos. Also, 
gluino will decay dominant into the stop and top quarks, which have long decay chains as 
well. Therefore, the LHC supersymmetry search constraints can be relaxed. The detailed 
LHC supersymmetry search constraints will be studied elsewhere. Moreover, because the 
LSP neutralino has relatively large Higgsino components due to small effective fi term, 
the XENON100 experiment gives strong constraint, for example, the spin-independent LSP 
neutralino-nucleon cross sections are larger than the XENON100 experiment upper bound 
in the bencharmark points IIIx^ in , IIIA, and IIIC [23J if i?-parity is conserved. This is 
another reason why we get better points in Scenario III than Scenario II. In Scenario III, 
the constraints from the LHC supersymmetry searhes and XENON100 experiment can be 
escaped, and the i?-parity violating \j k and \' i j k terms can increase (g M — 2)/2 and generate 
the neutrino masses and mixings. Therefore, Scenario III with i?-parity violation is more 
natural and realistic than Scenarios I and II. 



V. CONCLUSION 



We pointed out that as a solution to the SM fine-tuning problems, supersymmetry needs 
not to provide the dark matter candidate, i.e., i?-parity can be violated. Because the 
NMSSM can explain the Higgs boson mass and decay rates better than the MSSM, we 
considered three kinds of the NMSSM scenarios. To satisfy the phenomenological constraints 
and fit the experimental data, we studied the x 2 analyses for all three kinds of Scenarios, 
and showed that the Higgs boson mass and decay rates can indeed be explained very well. 
For the small x 2 values and fine-tuning around 2-3.7%, we obtained the viable parameter 
space with light (e.g. less than about 900 GeV) supersymmetric particle spectra only in 
Scenario III. With the small x 2 values and fine-tuning around 1-2%, we got the viable 
parameter space with relatively heavy (e.g. less than about 1.2 TeV) supersymmetric particle 
spectra. Especially, the best benchmark point is IIIA, which has almost minimal x 2 an d 
3.7% fine-tuning. The correct dark matter density can be realized in Scenario I as well. The 
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FIG. 2: The fitting results for Scenario I with relic density included in the x 2 ■ The red stars show 
the best-fitted benchmark point with minimal Xmin = 21.16. The green, blue, and black regions 
are respectively one, two, and three standard deviation regions with x 2 < Xmin + l^Ln + 4 an< ^ 
y 2 ■ +9 

generic features for the viable parameter space with smaller x 2 are that the light stop is 
around 500 GeV or smaller, the singlino and Higgsino are light chargino and neutralinos, 
the Wino-like chargino is heavy, and the Bino-like and Wino-like neutralinos are the second 
heaviest neutralino and heaviest neutralinos, respectively. Thus, we found that the LHC 
supersymmetry search constraints can be relaxed due to quite a few jets and/or leptons 
in the final states in Scenarios I and II, but the XENON100 experiment still gives strong 
constraint on the dark matter direct detection cross sections. Moreover, tan (3 is not large 
and the second lightest CP-even Higgs particle is SM-like so that the SM-like Higgs boson 
mass can be lifted. However, the extra contributions to the muon (g fJi — 2)/2 are smaller than 
three sigma low bound in general due to relatively small tan j3. With i?-parity violation, we 
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1.5 




FIG. 3: The fitting results for Scenario II with relic density smaller than the 95% C.L. upper 
limit. The red stars show the best-fitted benchmark point with minimal Xmin = 19.35. The 
green, blue, and black regions are respectively one, two, and three standard deviation regions with 

X < Xmin 1' Xmin ^' an( l X m in ^- 

can evade the LHC supersymmetry search and XENON100 experiment constraints, and the 
.R-parity violation superpotential term(s) may increase (g u — 2)/2 and explain the neutrino 
masses and mixings. Therefore, Scenario III with i?-parity violation is more natural and 
realistic than Scenarios I and II. 
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TABLE IV: Particle spectra (in GeV) and parameters for benchmark points in Scenario I. 
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75.7 


59.6 


e//i 


680 


377 


710 


489 


o 77 
^ V or 


1.08 


1.68 


1.53 


1.24 


V r 


679 


377 


709 


489 




1.34 


1.45 


1.42 


1.42 




376 


122 


510 


191 


R WW 


0.89 


1.09 


1.10 


1.09 




683 


383 


713 


492 


R 2 ZZ 


0.89 


1.09 


1.10 


1.09 


n 


372 


118 


509 


190 


R 2 Vbb 


0.01 


0.57 


0.72 


0.65 


r 2 


683 


383 


712 


492 


R 2 bb 


0.01 


0.49 


0.67 


0.74 


x 2 


19.35 


24.19 


23.86 


23.70 


R 2 TT 


0.00 


0.49 


0.66 


0.73 



TABLE V: Particle spectra (in GeV) and parameters for benchmark points in Scenario II. 
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Point 


Illy 2 


IIIA 


IIIB 


IIIC 


Point 


Illy 2 


IIIA 


IIIB 


IIIC 


M 


352 


431 


344 


337 


h 


365 


252 


184 


391 




489 


326 


500 


441 


h 


850 


647 


832 


795 


tan/3 


2.341 


1.853 


2.731 


2.039 


t 


925 


694 


852 


886 


A 


0.604 


0.589 


0.623 


0.614 


h 


814 


610 


793 


760 


K 


0.295 


0.364 


0.288 


0.318 




1012 


778 


1020 


929 


A 


-1063 


-372 


-1322 


-620 


ur/cr 


1078 


827 


1097 


985 


A x 


-329 


-161 


-267 


-9.93xl0~ 6 


ul/cl 


1083 


823 


1101 


990 


A K 


-128.134 


-628.214 


-0.877 


-1.459 


Ar/'sr 


1039 


798 


1056 


952 




144 


162 


150 


166 


cIl/sl 


1085 


825 


1103 


992 


Mi 


207 


136 


212 


186 


1 


120.5 


114.3 


119.6 


116.6 


M 2 


384 


255 


393 


346 


A 


126.9 


124.6 


124.9 


125.7 


M 3 


1090 


745 


1114 


990 


m 

o 


367.8 


342.7 


436.6 


382.5 


X? 


91 


86 


92 


107 


Ai 


154.3 


300.6 


157.1 


247.6 


xl 


-177 


160 


-184 


195 


A 2 


360.0 


335.1 


430.1 


376.3 




187 


-189 


191 


-197 




352.6 


330.2 


421.7 


369.3 




222 


232 


228 


221 




X 


X 


X 


X 


xl 


417 


310 


425 


384 


"-/J, L J 


1.893 


1.587 


2.207 


1.848 


xf 


126 


122 


133 


141 


a si (p) [10- 10 pb] 


X 


X 


X 


X 


xf 


416 


306 


425 


383 


Br( 6 ^) [10- 4 1 


3.586 


3.413 


2.560 


3.659 


g 


1138 


796 


1162 


1035 


Aft 


59.5 


27.0 


68.4 


44.1 


e/ fj, 


513 


502 


514 


475 


^ Vrir 


0.88 


0.93 


1.60 


1.16 


U T 


512 


502 


514 


475 




1.34 


1.46 


1.41 


1.42 




291 


395 


273 


289 


R WW 


0.81 


0.95 


1.08 


1.10 




517 


505 


519 


479 


R 2 ZZ 


0.81 


0.95 


1.08 


1.10 


n 


289 


395 


270 


288 


R 2 Vbb 


0.01 


0.10 


0.55 


0.37 


T2 


517 


505 


518 


479 


R2 bh 


0.02 


0.16 


0.48 


0.46 


x 2 


19.67 


21.31 


23.85 


20.53 


R 2 TT 


0.00 


0.14 


0.47 


0.44 



TABLE VI: Particle spectra (in GeV) and parameters for benchmark points in Scenario III. 
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